Light induced crystallization of cold atoms in a thin ID optical tube 
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Collective off resonant scattering of coherent light by a cold gas induces long-range interactions 
via interference of light scattered by different particles. In a ID configuration these interactions 
grow particularly strong for particles trapped along an optical nanofiber. We show that there exists 
a threshold pump laser intensity, above which the gas can be found in a crystalline, selfsustained 
order. In the nonabsorbing regime we determine the critical condition for the onset of order as 
well as the forms of particle density and scattered field patterns along the fiber above threshold. 
Surprisingly, there can coexist multiple stationary solutions with distinct density and field profiles. 

PACS numbers: 37.30.+i, 37.10.-x, 51.10.-|-y 



The astonishing experimental accomplishments in the 
optical control and manipulation of cold atomic gases in 
the past decade now allow to deterministically load them 
in extremely well controlled optical traps of almost any 
shape down to effectively zero temperature. A particu- 
lar fruitful example are periodic optical lattices in which 
many intriguing phenomena of solid state physics can be 
studied with unprecedented control[l]. 

In contrast to conventional solids, however, the spatial 
order and lattice geometry is fixed by the external lasers 
and does not appear from a self consistent dynamics of 
particle interactions. As the back-action of the atoms 
onto the confining light fields is generally negligible [2], lo- 
cal lattice perturbations do not propagate and long range 
interactions, which appear as phonons in solids are ab- 
sent. This changes for spatially tightly confined fields as 
in small optical resonators or optical micro-structures [3]. 

Experimentally it is very challenging to implement and 
load microtraps close to optical microstructures, where 
such atom field coupling is enhanced [4-6]. In an im- 
portant step Rauschenbeutel and coworkers, however, 
recently managed to trap atoms in an array of optical 
dipole traps generated by two color evanescent light fields 
alongside a tapered optical fiber [7] where the backaction 
of even a single atom on the propagating fiber field is 
surprisingly strong[3, 8]. This setup was improved with 
higher control and coupling by other groups recently [9]. 
With the atoms firmly trapped within the evanescent 
field of the fiber modes, field mediated atom-atom in- 
teraction and collective coupling to the light modes play 
a decisive role in this a setup[10]. 

Already a decade ago it was theoretically predicted[ll] 
and experimentally confirmed[12~14] that light scattering 
within optical resonators induces self-ordering of atoms 
in regular patterns maximizing collective coupling to the 
cavity mode [15, 16]. This transition can be directly mon- 
itored from the super-radiant light scattering [12, 14] and 
appears also at zero temperature in the quantum regime 
as phase transition from a superfluid to a supersolid[13]. 
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In an recent proposal Chang and coworkers predicted, 
that nanofiber-mediated infinite range dipole-dipole cou- 
pling can also induce stable regular patterns of laser il- 
luminated atoms trapped along the fiber[17]. The stable 
configurations, characterized by minimal dipole interac- 
tion energy, assume surprising configurations and exhibit 
characteristic collective light scattering. 

Here we develop a generalized many particle model to 
study the properties of such a light scattering induced 
crystallization of an laser illuminated ultracold gas in an 
elongated ID trap at finite temperature as schematically 
drawn in Fig.l. The atoms trapped parallel to the fiber 
are illuminated at right angle by a Gaussian laser beam 
of frequency w/,, sufficiently detuned from any atomic 
resonance so that spontaneous emission plays only a mi- 
nor role and their polarizability a has only a negligible 
imaginary part. 
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FIG. 1; Cigar-shaped atomic gas alongside optical nanoguide 

The laser field = (fL(x)e^"^^* -f c.c.)eL gives rise 
to the field Es scattered by the atoms. Neglecting po- 
larization effects and retardation, it can be written as 
Es = (fs(x, t)e~*'^^* -I- c.c.)ei, its envelope satisfying 
Helmholtz's equation 

V2£:, + fci(n|, + x)f. --fcix^L. (la) 

Here, ni?(x) denotes the refractive index profile of the 
optical fiber and x(^i^) — '^P^^'^')I^q is the suscep- 
tibility of the particles. The particle density distribu- 
tion p is obtained from the momentum integral /9(x, t) — 
J F{x,p,t)d'^p over the one-body distribution function 
F(x, p, t) proper to the gas. In the classical mean-field 
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limit it satisfies Vlasov's equation 

dF dF d OF _^ 

dt m 9x 9x 9p 

Here, Ut denotes the elongated dipole trap potential and 



$d = -a \£s + 



(2) 



the optical potential due to pump laser and fiber field. 
For strong radial confinement of the gas the one-body 
distribution approximately factorizes into a longitudinal 
and a transverse part, F(x, p, t) ~ f)Fj_(xj_, pj_), 

where F±_ is the Maxwell-Boltzmann distribution 
for the transverse degrees of freedom: _Fl := 

j3 denotes the inverse 



Z I ^ cxp 



thermal energy and U±_ the radial confining potential. 

In the absence of atoms the fiber supports only a sin- 
gle relevant TE mode, which - within the framework 
of scalar theory - is supposed to possess a propagation 
constant j3m and a normalized transverse mode function 
u{x, y) extending outside the fiber [7, 17]. The dominant 
forces on the particles along z are due to photon scatter- 
ing into and out of the fiber. As long as the total atomic 
susceptibility stays small even for large particle numbers, 
the radial fiber mode function is only weakly perturbed 
and we can set iPs(x, t) ~ \fAE{z, t)u{x, y) with the cross 
section ^4 := (/ (Px±(Pp±u^F±) . Integrating over the 
transverse degrees of freedom we finally arrive at an ef- 
fective description of the longitudinal dynamics 



{Pi + klx) E = -klxEL 



(3a) 



dt m dz 



dz 



{U-a[\E\^ + 2ELEr]) 



dl 
dpz 



0, 

(3b) 

where Er is the real part of E. The effective laser field 
is given hy E^ := V A J d'^x±d'^pj^£i^uFj_ and the local 
atom-fiber coupling is governed by the effective suscepti- 
bility 



X{z,t) := 



f 

J — c 



f{z,Pz,t)dpz 



(4) 



proportional to the atomic line density. To ensure phys- 
ically consistent solutions, equation (3a) for the electric 
field has to fulfill Sommerfeld's radiation conditions 



dE 

dz 



±iPrnE{z, t), Z ^ ±00, 



(5) 



which ensure outgoing waves at infinity. 

Equations (3) can at least numerically be solved di- 
rectly. At this point we nevertheless restrict ourselves 
to such stationary solutions, which correspond to a gas 
in thermal equilibrium [18], where the particles are dis- 
tributed according to the Maxwell-Boltzmann distribu- 
tion 



f{z,Pz) 



-l3(pl/2m+U) Ba{\E\^+2EpE,) 



(6) 



Substituting the effective susceptibility (4) obtained 
from a thermal distribution (6) into the effective 
Helmholtz equation (3a) leads to a highly nonlinear equa- 
tion, which determines the selfconsistent electric field. 
We assume here that the external part of the poten- 
tial, f/(z), can be harmonically approximated by U{z) = 
^muj'^z'^ within the (very large) thermal extension of the 
gas Iz := 2kBT/muj1 ^ /3~^, and we will ignore the de- 
pendence of the driving laser amplitude on position, i.e. 
El{z) = El. 

With these assumptions one sees that the zero field 
solution E{z) w and 



fiz,P^) = fo{z,Pz) = Z 



(7) 



always solves Helmholtz's equation, which we will call 
the normal phase. 

Before we investigate the possibility of selfordered so- 
lutions, let us examine the dynamical stability of the 
normal phase, using the Helmholtz- Vlasov equations (3), 
linearized around the unordered state 



— +0^E--El^ 



h{z,Pz,t)dpz, (8a) 
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^ Pzdfi 
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dUdfi 
dz dpz 



= -2a 



dELEr dfo 

dz dpz 



(8b) 



Here, f\{z,Pz,t) represents a small deviation from the 

normal phase, such that the complete distribution func- 
tion is given by f{z,pz,t) = fo{z,Pz) + fi{z,Pz,t). We 
seek solutions to (8), which are of the form 



fi{z,Pz,t) = ip{z,pz)e'* + 4}*{z,pz)e- 



*t 



E{z,t) = a{z)e'* + b*{z)e 



(9a) 



(9b) 



In this Ansatz, (tjj,a,b) will be referred to as normal 
mode and s = 'y + iui as the complex valued mode pa- 
rameter. Whenever 7 > 0, the deviation from equilib- 
rium defined by the normal mode increases exponentially 
in time, causing the eventual destruction of the normal 
phase. The existence of such a mode implies dynamical 
instability. With the help of asymptotic expansions, we 
find that if (9) is to solve (8) and satisfy the radiation 
conditions, the mode parameter must satisfy Dn{s) = 
for some n G Z, where 



2 ^,2 r r°° p2 



Z)„(s) = (2n + l)7r-hi^ J J 



Ejdfo/dpz 

S + iPmVz 



dpzdz. 
(10) 

This condition first demands w = 0. Introducing 
the collective coupling parameter ^0 and effective pump 
strength e 



/- kp 

Co •- 
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(11) 
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FIG. 2: Properties of three coexisting selfordered solutions on 
different branches n = ((lo) — (Ic)), n = 1 ((2a) — (2c)) and 
n = 2 ((3a) — (3c)) in the strong collective coupling regime 
for C,o — ISO/tt and e = 9.5£c- The upper plots (a) show 
the atomic density along the fiber, while the plots labeled 
by (c) depict the local fraction of right (blue curve) and 
left A'^- (red curve) traveling photons, the number of zeros 
being equal to the branch number n. The plots labeled (b) 
depict the envelopes of the two parts of the optical potential. 
Blue indicates the potential from scattering between pump 
laser and fiber, the brown area shows the part due to photon 
redistribution within the fiber. For n — the latter dominates 
in the center of the cloud causing a density modulation with 
a period of one half of Am = 2ti/ Pm- This region is broken up 
into three for n = 1 and into five, for n = 2. 



we then see that there exists at least one normal mode 
with a positive growth rate 7 > 0, if and only if 

^>^^-- (12) 

Hence as our first important result we see that be- 
yond this pump threshold (12) the normal phase ceases 
to be stable and no longer represents a physical state of 
the system. We therefore study other solutions of (3) 
for e > Ec- Fortunately, we can obtain approximate so- 
lutions with the help of perturbation theory in the weak 
collective coupling regime, Co ^ 1; as well as in the strong 
collective coupling regime Co 3> 1. The real and imagi- 
nary parts of the electric field envelope Er,Ei 

Er,i{z) = ar^i{z) C0S(^„Z -I- (j>r^i{z))EL, (13) 

can be assumed to behave almost harmonically with 
phases (t>n4'i a-nd amplitudes ar,ai that vary slowly on 
a scale defined by This amplitudes and phases 

are directly related to the complex amplitudes of the 
more familiar decomposition into left and right run- 
ning waves, E = (£^+6*'^™^ + E_e-'P^'')EL via E± = 
Ore^"^'' -I- iflie**'^' . A perturbative analysis of the steady 
state Helmholtz equation reveals, that Q := ''^ ' = 
\ ( I E+ p -I- 1 p ) , proportional to the sum of locally right 
going and left going photons, is spatially constant and 
will thus serve us as order parameter. 

Let us first consider the weak collective coupling 
regime, where scattering from the laser into the fiber 
and vice versa dominates over scattering within the fiber. 
The demand to satisfy the radiation conditions (5) leads 
to the equation determining this order parameter 

2Co/i (2e%/e) = (1 + 2n)\/e/o (2e%/e) , n € No, 

(14) 

where Ik denotes the fcth modified Bessel function of the 
first kind. 

As quite a surprise one observes that the solution is 
not unique and for Ec < £ < j there exist n different 
solutions Qn as shown in the example of figure 3. As the 
phase of the outgoing scattered light is not determined, 
each solution of (14) corresponds to an infinite family 
of solutions of (3) with the periodic density modulation 
slightly displaced under a slow envelope. The appearance 
of this degeneracy of course corresponds to breaking of 
a continuous symmetry. It is important to note that the 
point of first appearance of an ordered family of solu- 
tions, s — Ec, exactly coincides with the point where the 
normal phase becomes unstable. All this is confirmed 
by a numerical solution of the underlying equations (3) . 
Focusing on the region close to the first branching point, 
we find the behavior of the order parameter to be given 

by 

e-£-£„ Co<i (15) 
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FIG. 3: Branches of the order parameter vs energy ratio ac- 
cording to (14) and (18) for n = 0, . . . , 3 in the weak collective 
couphng regime (a), with (^o = 0.05 and in the strong collec- 
tive coupling regime (b), with = TS/vr. Note that in the 
latter regime all branches converge and there may be even 
more than n families of solutions for Ec < e < (1 -I- 2n)ec. 



and 



e - (e - Ec) 



1/2 



Co 



1. 



(16) 



Curiously, if > Ij the order parameter is discontinuous, 
exhibiting a jump of magnitude 



Ae - 16(Co - 1). 



(17) 



The results of the perturbation theory hold only as long 
as ^ 1 and the predicted jump soon violates this as- 
sumption. Let us therefore turn to the strong collective 
coupling regime. Again, each possible value of the order 
parameter corresponds to an infinite family of stationary 
solutions and any nonzero value for Q satisfies 



where P{&) :~ 4+eee ^^'^ ^ denotes the complete el- 
liptic integral of the first kind. In this regime, there 
are at least n families of ordered solutions whenever 
1 < 2Co£ < 1 + 2rt. For an example see figure 3. Further- 
more, we find that according to perturbation theory, the 
order parameter is bounded by O < 4. Figure 2 depicts 
some properties of the solutions corresponding to possi- 
ble values of the order parameter. It should be noted 
that while we have here shown the existence of multiple 
families of selfordered thermal solutions of the Vlasov- 
Helmholtz equations, nothing has been said about their 
dynamical stability properties, which remains an open 
problem. 

Certain phenomena called optical binding, similar to 
the ordering discussed in this work, have been observed 
with laser illuminated small beads in liquids[19] in ID 
and even 2D geometries as well. More interestingly, in 
analogy with cavity induced selfordering we expect a cor- 
responding phase transition also in the zero temperature 
limit, where the threshold is determined by the recoil 
energy and effective particle interaction strength replac- 
ing thermal energy. This should have distinctly different 
properties as compared to a ID prescribed optical lattice. 
Obviously, even without the presence of the fiber, the in- 
terference of the scattered field by two distant atoms can 
induced long range forces and ordering. As part of the 
scattered field is lost, this leads to a higher threshold, but 
to some extend the particles itself can guide the scattered 
light. In the limit of sufficient initial particle density the 
trapped atoms themselves could form a fiber like guide 
for the scattered light enhancing long range interactions. 
One could speculate that the trapping field for the ID 
confinement could be unnecessary and be provided by 
the light guided by the atoms. The predicted ordering 
phenomena could also be related to recent observations of 
collective scattering in dense atomic vapors [20, 21] where 
light induced selfordering mechanisms play an important 
role. 
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Appendix 

In this appendix wc indicate, how the results of the 
present work can be obtained in a systematic manner by 
means of canonical perturbation theory, not restricting 
ourselves to thermal equilibrium. Helmholtz's equation 
with a stationary gas density p in the source term reads 



El)P- 



(19) 



It is important to note that for any stationary solu- 
tion of Vlasov's equation (3b), which necessarily takes 
the form f{z,pz) — F{p'^/2m + + U), where $d = 
— a(|£^|^ + 2ELEr), we have for the corresponding den- 
sity p{z) = g[^d + U] for some real function g. Let G 
denote the antiderivative of g, i.e. G' = g, and let us 
define the canonical momenta Ilr^i := d^E^^i. Then, the 
above equations for the real and imaginary parts of the 
selfconsistent field are the canonical equations of the fol- 
lowing Hamiltonian 



H 



ILl+lfi+/3l{E',+Ef) 



u] 



(20) 

which has the physical interpretation of being propor- 
tional to the energy density of the scattered electric field. 
It can be written as the Hamiltonian of two uncoupled 
resonant oscillators, Hq, plus a coupling, H = Hq + Hi, 
where 



2eoA 



G[^d + U], 



(21) 



which part we shall consider a perturbation. Let us 
also introduce action-angle variables {Jr, Jijipr^i^i) for 
the unperturbed part via the prescription 



Er,i = W -7r^sin('^r,i) 



/3„ 



(22a) 



(22b) 



Then we find the transformed Hamiltonian as H = 

Ho{J) + eHi{J,ip,ez), where now Hf){J) — f3m{Jr + Ji) 
and we have introduced a purely formal expansion pa- 
rameter e. The explicit dependence on z is due only to 
the external trap U{z) and thus weak. Next wc will seek 
a canonical transformation, effected by the mixed vari- 
ables generating function S = S{J, ip, z) to new variables 
{J,'ip) = {Jr,Ji,'il>r,i'i) according to 



(23) 



and leading to the new Hamiltonian K = H + which 

we wish to make as simple as possible. To exploit the 
assumed smallness of Hi , we express the generating func- 
tion as 

S = J.^ + eSi+e^S2 + .... (24) 
Likewise expanding K — Kq + eKi + . . . , we find 

Kq{J) = Ho{J), (25a) 



Ki{J,ij,€z) = V -^+Hi{J,^,ez), (25b) 

and so forth. Setting e — 1 and introducing the unper- 
turbed trajectories 'ipr,i{C) '■= '4'r,i + C) we find 



Si{J,^ij{0,z)-SiiJ,^P,z) 
rC 

A" ' 



^ [Ki{J,^{C),z) - Hi{J,^{C),z)] dC. (26) 

'm Jo 



Imposing the periodicity of the generating function in the 
angles leads to the solvability condition 



r [Ki{J,i;{C),z) - HiiJ,i;{0,z)] dC = 0, (27) 
Jo 

which must be satisfied by any potential Ki . Now we use 
our freedom to choose Ki to be a function of the phase 
difference only, Ki = Ki{J, A, z), where A := -ipr — 't/Ji. 
Then it is easy to see, that (27) demands 

Ki{J, A, z) := ^ j^^ Hi{J, C + A, C, z)dC. (28) 

In this way we end up with the new and simplified 
Hamiltonian, to first order in the perturbation given by 
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K{J,ijj,z) = Ho{J) + Ki{J,A,z). As this new Hamil- 
tonian depends, by construction, only on the phase dif- 
ference A, we immediately obtain an integral of motion, 
8 := Jr + Ji, proportional to the order parameter, which 
is the total number of quanta of excitation in the left and 
right propagating field components. Therefore, the num- 
ber of degrees of freedom is reduced to just two, namely 
A and D := Jr — Ji, whose equations of motion are the 



canonical equations pertaining to the Hamiltonian 

H{D,A,z):=2k^{^^,^^,A,z^ (29) 

and can be solved in the weak and strong collective cou- 
pling limits respectively. 



